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<^ ■ Abstract 

The diffraction of stochastic point sets, both Bernoulli and Markov, and 
| of random tilings with crystallographic symmetries is investigated in rigorous 

terms. In particular, we derive the diffraction spectrum of ID random tilings, 
■ of stochastic product tilings built from cuboids, and of planar random tilings 

based on solvable dimer models, augmented by a brief outline of the diffraction 



ON , sure theoretic approach to mathematical diffraction theory which underlies the 

unique decomposition of the diffraction spectrum into its pure point, singular 
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^ Introduction 



from the classical 2D Ising lattice gas. We also give a summary of the mea- 
sure theoretic approach to mathematical cliff 
unique decomposition of the diffraction spec 
continuous and absolutely continuous parts. 
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The diffraction theory of crystals is a subject with a long history, and one can safely 
say that it is well understood [^, |T^]. Even though the advent of quasicrystals, with 
their sharp diffraction images with perfect non-crystallographic symmetry, seemed to 
question the general understanding, the diffraction theory of perfect quasicrystals, in 
terms of the cut and project method, is also rather well understood by now, see |26], |27j 
and references therein. It should be noted though that this extension was by no means 
automatic, and required a good deal of mathematics to clear up the thicket. More 
recently, this has found a general extension to the setting of locally compact Abelian 
groups f)!| which can be seen as a natural frame for mathematical diffraction 
theory and covers quite a number of interesting new cases f|. 
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Another area with a wealth of knowledge is the diffraction theory of imperfect 
crystals and amorphous bodies |2(| ^8|, but the state of affairs here is a lot less 
rigorous, and many results and features seem to be more or less folklore. For example, 
the diffraction of simple stochastic systems, as soon as they are not bound to a lattice, 
is only in its infancy, see [|J for some recent addition to its rigorous treatment. This 
does not mean that one would not know what to expect. However, one can often only 
find a qualitative argument in the literature, but no proof. May this be acceptable 
from a practical angle, it seems rather unsatisfactory from a more fundamental point 
of view. In other words, the answer to the question which distributions of matter 
diffract is a lot less known than one would like to believe, compare the discussion in 



2|, Sec. 6] and also |IJ |7j]. 

Note that this question contains several different aspects. On the one hand, one 
would like to know, in rigorous terms, under which circumstances the diffraction image 
is well defined in the sense that it has a unique infinite volume limit. This is certainly 
the case if one can refer to the ergodicity of the underlying distribution of scatterers 



Toj , |26| , p3| , in particular, if their positional arrangement is linearly repetitive p| . 
However, this is often difficult to assess in situations without underlying ergodicity 
properties, see |5[] for an example. On the other hand, even if the image is uniquely 
defined, one still wants to know whether it contains Bragg peaks or not, or if there is 
any diffuse scattering present in it. 

This situation certainly did not improve with the more detailed investigation of 
quasicrystals, e.g. their less perfect versions, and in particular with the study of the 
so-called random tilings [[16], |47j . Again, there is a good deal of folklore available, 
and a careful reasoning based upon scaling arguments (compare ^j) seems to 
give convincing and rather consistent results on their diffraction properties. However, 
various details, and in particular the exact nature of the diffraction spectrum, have 
always been the topic of ongoing discussion, so that a more rigorous treatment is 
desirable. It is the aim of this contribution to go one step into this direction, and to 
extend the analysis of [[| on generalized lattice gases to the case of certain Markov 
type systems as they appear in the theory of random tilings. We will not be able 
to answer the real questions concerning those tilings relevant for quasicrystals, but 
we still think that the results derived below are a worth-while first step. Even this 
requires a bunch of methods and results which are scattered over rather different 
branches of mathematics and mathematical physics. It is thus also one of our aims to 
recollect the essential aspects and references, tailored for what we need here and for 
future work in this direction. 

Let us summarize how the article is organized. We start with a recapitulation 
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of the measure theoretic setup needed for mathematical diffraction theory, where we 



essentially follow Hof |26], |27|, but adapt and extend it to our needs. We will be a little 
bit more explicit here than needed for an audience with background in mathematics or 
mathematical physics, because we hope that the article becomes more self-contained 
that way, and hence more readable for physicists and crystallographers who usually 
do not approach problems of diffraction theory in these more rigorous terms. We 
consider this as part of an attempt to penetrate the communication barrier. We then 
investigate several ID systems, notably Markov systems and ID random tilings, and 
derive their spectral properties. This is followed by an intermediate discussion of 
stochastic product tilings in arbitrary dimensions which already indicates that the 
appearance of mixed spectra with pure point, singular continuous and absolutely 
continuous parts is generic, though it also shows that its meaning in more than one 
dimension will have many facets. 

The next Section then deals with the main results of this article, the derivation 
of the diffraction spectrum of certain crystallographic random tilings in the plane, 
namely the domino and the rhombus (or lozenge) tiling. This requires some adapta- 
tion of results from the theory of Gibbs states to special hard-core lattice systems. 
Again, we explain that in slightly more detail than necessary from the point of view 
of mathematical physics in order to enhance self-containedness and readability. In 
both cases, the basic input of the explicit result has long been known in statistical 
mechanics, but the interest in the diffraction issue is rather recent. We also briefly 
comment on the diffraction of an interactive lattice gas based upon the classical 2D 
Ising model and its implications. The discussion addresses some open questions and 
what one should try to achieve next. 

Recollections from mathematical diffraction theory 

Diffraction problems have many facets, but one important question certainly is which 
distributions of atoms lead to well-defined diffraction images, and if so, to what kind 
of images. This is a difficult problem, far from being solved. So, one often starts, 
as we will also do here, by looking at "diffraction at infinity" from single-scattering 
where it essentially reduces to questions of Fourier analysis Sec. 6]. This is also 
called kinematic diffraction in the Fraunhofer picture fT2|, and we are looking into 



the more mathematical aspects of that now. Mathematical diffraction theory, in turn, 
is concerned with spectral properties of the Fourier transform of the autocorrelation 
measure of unbounded complex measures. Let us therefore first introduce and discuss 



the notions involved. Here, we start from the presentation in [26, 27 where the linear 
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functional approach to measures is taken, compare [I3| for details and background 
material. We also introduce our notation this way. 

Let K be the space of complex- valued continuous functions with compact support. 
A (complex) measure fi on M n is a linear functional on K, with the extra condition 
that for every compact subset K of M. n there is a constant a K such that 

Hf)\ < a K \\f\\ (1) 

for all / G K, with support in K; here, ||/|| = sup^^- \f{x)\ is the supremum norm of 
/. If p is a measure, the conjugate of \i is defined by the mapping / — > p(f). It is 
again a measure and denoted by p. A measure /x is called real (or signed), if p = p, 
or, equivalently, if is real for all real-valued / G K. A measure p is called positive 
if //(/) > for all / > 0. For every measure /i, there is a smallest positive measure, 
denoted by such that |/x(/)| < |/i|(/) for all non-negative / G /C, and this is called 
the absolute value (or the total variation) of //. 

A measure \i is bounded if |/x|(R n ) is finite (with obvious meaning, see below), other- 
wise it is called unbounded. Note that a measure \i is continuous on K with respect to 
the topology induced by the norm if and only if it is bounded [0, Ch. XIII. 20]. In 
view of this, the vector space of measures on M n , A4{R n ), is given the vague topology, 
i.e. a sequence of measures {fi n } converges vaguely to p if lim n _^ ooj u n (/) = p(f) in C 
for all / G K. This is just the weak-* topology on A4(lR n ), in which all the "stan- 
dard" linear operations on measures are continuous, compare |46| , p. 114] for some 
consequences of this. The measures defined this way are, by proper decomposition 



14] , Ch. XIII. 2 and Ch. XIII. 3] and an application of the Riesz- Markov representation 



theorem, see 0, Thm. IV. 18] or || Thm. 69.1], in one-to-one correspondence with 
the regular Borel measures on M n , wherefore we identify them. In particular, we write 
fi(A) (measure of a set) and /i(/) (measure of a function) for simplicity. 

For any function /, define / by f(x) := f(—x). This is properly extended to 

measures via p(f) '■= fi>(f)- Recall that the convolution /i * v of two measures /i and 
v is given by \i * v(f) := f f(x + y)p{dx)v{dy) which is well-defined if at least one 
of the two measures has compact support. For R > 0, let B R denote the closed ball 
of radius R with centre 0, and vo\(B R ) its volume. The characteristic function of a 
subset A C M. n is denoted by 1 A . Let p, R be the restriction of a measure \i to the ball 
B R . Since fi R then has compact support, 

7 * := voW^*^ (2) 

is well defined. Every vague point of accumulation of j R , as R — » oo, is called an 
autocorrelation of //, and as such it is, by definition, a measure. If only one point 
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of accumulation exists, the autocorrelation is unique, and it is called the natural 
autocorrelation. It will be denoted by 7 or by 7^ to stress the dependence on \i. One 
way to establish the existence of the limit is through the pointwise ergodic theorem, 
compare |15|], if such methods apply. If not, explicit convergence proofs will be needed, 



as is apparent from known examples [Q and counterexamples |40[ . 



Note that Hof |2(J uses cubes rather than balls in his definition of j R . This simpli- 
fies some of his proofs technically, but they also work for balls which are more natural 
objects in a physical context. This is actually not important for our purposes here. 
One should keep in mind, however, that the autocorrelation will, in general, depend 
on the shape of the volume over which the average is taken — with obvious meaning 
for the experimental situation where the shape corresponds to the aperture. To get 
rid of this problem, one often restricts the class of models to be considered and defines 



the limits over van Hove patches, thus demanding a stricter version of uniqueness |50 
Sec. 2.1]. 

The space of complex measures is much too general for our aims, and we have to 
restrict ourselves to a natural class of objects now. A measure fi is called translation 
bounded |I| if for every compact set K C lR n there is a constant b K such that 

sup \/j,\(K + x) < b K . (3) 

xe«. n 

For example, if A is a point set of finite local complexity, i.e. if the set A = A — A of 
difference vectors is discrete and closed, the weighted Dirac comb 

:= ^w(x)5 X) (4) 

where S x is Dirac's measure at point x, is certainly translation bounded if the w(x) are 
complex numbers with sup xgA |u>(a;)| < 00. This is so because A discrete and closed 
implies that G A is isolated and the points of A are separated by a minimal distance, 
hence A is uniformly discrete. Translation bounded measures \x have the property that 
all 7 R are uniformly translation bounded, and if the natural autocorrelation exists, it 
is clearly also translation bounded |26], Prop. 2.2]. This is a very important property, 
upon which a fair bit of our later analysis rests. Note that such a restriction is neither 
necessary, nor even desirable (it would exclude the treatment of gases and liquids), 
but it is fulfilled in all our examples and puts us into a good setting in all cases where 
we cannot directly refer to pointwise ergodic theorems. Let us finally mention that 
different measures can lead to the same natural autocorrelation, namely if one adds 
to a given measure \i a sufficiently "meager" measure u, see ]26|, Prop. 2.3] for details. 
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In particular, adding or removing finitely many points from A, or points of density 0, 
does not change 7, if it exists. 

Let us focus on the Dirac comb uo = u A from (f|), with A of finite local complexity, 
and let us assume for the moment that its natural autocorrelation 7^ exists and is 
unique (Lagarias and Pleasants construct an example where this is not the case ffil). 



A short calculation shows that uj k = J2xeA w ( x )d-x- Since 5 X * 5 y = 5 x+y , we get 

7. = , (5) 

where the autocorrelation coefficient v{z), for z G A, is given by the limit 

u(z) = lim 1 Yl w{y)w(z-y), (6) 

z-yeA 

where An —An Br. Conversely, if these limits exist for all z G A, the natural 
autocorrelation exists, too, because A is discrete and closed by assumption, and (|) 
thus uniquely defines a translation bounded measure of positive type. This is one 
advantage of using sets of finite local complexity. 

We now have to turn our attention to the Fourier transform of unbounded measures 
on R™ which ties the previous together with the theory of tempered distributions [[54 
see ]T|, [53| for extensions to other locally compact Abelian groups. 



Let <S(R n ) be the space of rapidly decreasing functions |54|, Ch. VII. 3], also called 
Schwartz functions. By the Fourier transform of a Schwartz function G 5(M n ) we 
mean 

{T(j)){k) = 4>(k) := / e- 2 *** 4>(x)dx (7) 



which is again a Schwartz function [p4| , [4q| . Here, k ■ x is the Euclidean inner product 
of M. n . The inverse operation is given by 

ij(x) = [ e 2mx ' k i)(k)dk . (8) 

The Fourier transform T is thus a linear bijection from 5(lR n ) onto itself, and is 



bicontinuous j|6], Thm. IX. 1]. Our definition (with the factor 2tt in the exponent) 
results in the usual properties, such as <fi = <fi and tp = ip. The convolution theorem 
takes the simple form 0i * 02 = 0i ' 02 where convolution is defined by 

0i* 02 (x) := / 0i(x - y)(f) 2 {y)dy . (9) 
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Let us also mention that T has a unique extension to the Hilbert space L 2 (R n ), often 
called the Fourier-Plancherel transform, which turns out to be a unitary operator of 
fourth order, i.e. JF 4 = Id. This is so because (jF 2 0)(x) = 0(— x), see (|9| for details. 



Finally, the matching definition of the Fourier transform of a tempered distribution 
5| Te S'(R n ) is 

f (0) := T{4>) (10) 

for all Schwartz functions 0, as usual. The Fourier transform is then a linear bijection 
of iS^M" - ) onto itself which is the unique weakly continuous extension of the Fourier 



transform on 5(IR n ) f46[ , Thm. IX. 2]. This is important, because it means that weak 
convergence of a sequence of tempered distributions, T n — > T as n — > oo, implies weak 
convergence of their Fourier transforms, i.e. T n — > T . 

Let us give three examples here, which will reappear later. First, the Fourier 
transform of Dirac's measure at x is given by 

S x = e -2*f*-* (11) 

where the right hand side is actually the Radon-Nikodym density, and hence a function 
of the variable k, that represents the corresponding measure (we will not distinguish 
a measure from its density, if misunderstandings are unlikely). Second, consider the 
Dirac comb u r = XLer $e °f a lattice r C K" (i.e. a discrete subgroup of MJ 1 such 
that the factor group M n /r is compact). Then, one has 

ui v = dens(r) • u r , , (12) 

where dens(r) is the density of T, i.e. the number of lattice points per unit volume, 
and T* is the dual (or reciprocal) lattice, 

T* := {y G R n | x ■ y G Z for all x G T} . (13) 



This is Poisson's summation formula for distributions j54], p. 254] and will be central 
for the determination of the Bragg part of the diffraction spectrum. Finally, putting 
these two pieces together, we also get the formula 

^e~ 2 ^ = dens(r) • £ 5 y , (14) 

to be understood in the distribution sense. 

If a measure /i defines a tempered distribution T M by T M (0) = /z(0) for all G 
5(IR n ), the measure is called a tempered measure. A sufficient condition for a measure 
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to be tempered is that it increases only slowly, in the sense that f(l + \x\) e \p\(dx) < 



oo for some £ G N, see [54, Thm. VII. VII]. Consequently, every translation bounded 
measure is tempered - and such measures form the right class for our purposes. We 
will usually not distinguish between a measure and the corresponding distribution, 
i.e. we will write p, for T^. The Fourier transform of a tempered measure is a tempered 
distribution, but it need not be a measure. However, if p is of positive type (also called 
positive definite) in the sense that p(4> * (f>) > for all G 5(R ri ), then p is a positive 



measure by the Bochner-Schwartz Theorem [|6, Thm. IX. 10]. Every autocorrelation 
7 is, by construction, a measure of positive type, so that 7 is a positive measure. This 
explains why this is a natural approach to kinematic diffraction, because the observed 
intensity pattern is represented by a positive measure that tells us which amount of 
intensity is present in a given volume. 

Also, taking Lebesgue's measure as a reference, positive measures p admit a unique 
decomposition into three parts, 

p = fi m + p sc + p ac , (15) 

where pp, sc and ac stand for pure point, singular continuous and absolutely continu- 
ous, see [[K], Sec. 1.4] for background material. The set P = {x | /i({x}) 7^ 0} is called 
the set of pure points of /i, which supports the so-called Bragg part p of p. Note that 
P is at most a countable set. The rest, i.e. p — p , is the "continuous background" of 
/1, and this is the unambiguous and mathematically precise formulation of what such 
terms are supposed to mean. Depending on the context, one also writes 

V = Vpp + Vcont = Using + Vac , ( 16 ) 

where p cont = /i sc + fi ac = fi — fJ, pp is the continuous part of p, (see above) and p sing = 
p pp + p sc is the singular part, i.e. p sing (S) = for some set S whose complement has 
vanishing Lebesgue measure (in other words, p sing is concentrated to a set of vanishing 
Lebesgue measure). Finally, the absolutely continuous part, which is called diffuse 
scattering [^TJ in crystallography, can be represented by its Radon-Nikodym density 
f46| , Thm. 1.19] which is often very handy. Examples for the various spectral types can 



easily be constructed by different substitution systems, see |45| and references therein 
for details. Later on, we shall meet a simple example in the context of stochastic 
product tilings where all three spectral types are present, though their meaning will 
need a careful discussion. 

Hof discusses a number of properties of Fourier transforms of tempered measures 



2q , |27| . Important for us is the observation that temperedness of p together with 



positivity of p implies translation boundedness of p |2^, Prop. 3.3]. So, if p, is a 
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translation bounded measure whose natural autocorrelation 7^ exists, then 7^ is also 
translation bounded (see above), hence tempered, and thus the positive measure 7^ 
is also both translation bounded and tempered. This is the situation we shall meet 
throughout the article. 

In what follows, we shall restrict ourselves to the spectral analysis of measures /x 
that are concentrated on uniformly discrete point sets. They are seen as an idealization 
of pointlike scatterers at uniformly discrete positions, in the infinite volume limit. 
The rationale behind this is as follows. If one understands these cases well, one 
can always extend both to measures with extended local profiles (e.g. by convolution 
of uj with a smooth function of compact support or with a Schwartz function) and 
to measures that describe diffraction at positive temperatures (e.g. by using Hof's 



probabilistic treatment p8[). The treatment of gases or liquids might need some 
additional tools, but we focus on situation that stem from solids with long-range 
order and different types of disorder, because we feel that this is where the biggest 
gaps in our understanding are at present. 



Illustrative results in one dimension 

The simplest cases to be understood are those in one dimension. We start with some 
examples obtained from stationary stochastic processes and then derive in detail the 
diffraction properties of ID random tilings. The language and methods of this section 
closely follow those of classical ergodic theory because very good literature is available 
here @. 



Bernoulli and Markov systems 

Let us start with a Bernoulli system, i.e. with a lattice gas without interaction. 

Proposition 1 Consider the stochastic Dirac comb uj = Yl m ez r l{ rn )fim where r](m) 
is a family of i.i.d. random variables that can take any of the n complex numbers 
h x , . . . ,h n (assumed pairwise different), with attached probabilities p l7 . . . ,p n , p i > 0. 

Then, the autocorrelation 7^ of uj exists with probabilistic certainty^ and has the 
form 7^ = Xlmez u { m )^m with autocorrelation coefficients 

, ((\h\ 2 ), ifm = 0, , s 

vim) = I Xl 1 n J ' 17 

\\(h)\\ z/m^0, 



^ere and in the sequel, assertions of probabilistic certainty always refer to the invariant measure 
of the corresponding stochastic process. 
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where {\h\ 2 ) = Y^i=iPi\^i\ 2 an( ^ C 1 ) = SILi^iiV Consequently, the diffraction mea- 
sure is, with probability one, Z-periodic and given by 



iwEm + ((i^i 2 )-iwi 2 ) • 

v m£Z / 

Proof: This is a straight-forward application of Birkhoff 's pointwise ergodic theorem 



44] , Thm. 2.3], applied to the case of a Bernoulli system as described in |44j, Sec. 1.2 C]. 
One identifies the possible Dirac combs with the corresponding bi-infinite sequences 
x = (xj) i6 2 of the Bernoulli process given above. Then, v{m) is the orbit average of 
the function f(x) = x x m under the action of the shift which almost surely exists and 
equals the average of / over the invariant measure, because the process is ergodic. □ 

The diffraction thus consists of a pure point part (Bragg peaks) that is the one 
of the regular lattice Z multiplied by the absolute square of the average scattering 
strength and an absolutely continuous part (diffuse scattering) which is constant in 
this case (hence it is "white noise"), as one would expect for a Bernoulli process. Note 



that the entropy density of this ensemble is given by s = — X/iLi .Pi l°s(Pt)> see B 
Ch. 5.3, Ex. 3.4]. One could, alternatively, refer to the strong law of large numbers, 
under slightly different assumptions. This would then also give a generalization to 
higher dimensions, and to regular point sets beyond lattices, see j|] for a detailed 
account of this. 

Let us now turn our attention to stochastic systems with interaction. Let M be a 
Markov (or stochastic) matrix, i.e. M = (Mjj) 1<ij<n with My > and YTj=i Mij = 1. 
We assume that M is primitive, so some power of M has strictly positive entries 
only. As a consequence, the Perron-Frobenius (PF) eigenvalue Ai = 1 is unique, and 
all other eigenvalues Aj of M have absolute value |A$| < 1. The corresponding right 
eigenvector is (1,1,... ,1)*, while the left eigenvector p = pM, p = (p 1; . . . ,p n ), 
defines the stationary state. Due to the primitivity of M, we can choose p i > 
and statistical normalization, Y^i=xPi = 1- Let II = diag(p 1; . . . ,p n ) which is thus 
invert ible. 

In view of the applications we have in mind, we want to consider a Markov process 
that gives the same result for an operation in reverse direction. That is to say we 
restrict ourselves to reversible processes, i.e. to stochastic matrices M with 

ITM = MTI. (19) 

Since M is also primitive, all p i > and P = Il 1//2 is well defined and non-singular, 
P = diag( v /p^", . . . , y/p^)- Now, S = PMP~ X is a real symmetric matrix and can 
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thus be diagonalized by an orthogonal matrix. The eigenvalues of S and M coincide 
and are real; we denote them by Ai = 1, A 2 , . . . , A n , where |Aj| < 1 for all i > 2. If 
{bi}i<i< n is the corresponding orthonormal basis, the spectral decomposition of S is 

n 

S = l&iX&il© M-M&il = So ©Si, (20) 

i=2 

where we use Dirac's bra-ket notation for the standard Hermitian scalar product of C n . 
We embed MJ 1 into C n because we deal with complex scattering strengths later. Note 
that the decomposition on the right hand side of ( p0|) is into a projector (first term, 
So) and a contraction, i.e. \Six\ < \x\ for all x G C n . Also, we have S Si = SiS = 
and, in the standard basis of C n , So is explicitly given by S = (y^pj)i<ij< n - 

Proposition 2 Consider the stochastic Dirac comb uj = Y2m£Z r l{' m )$m where r](m) is 
a family of random variables that take values out of the n complex numbers h 1 , . . . ,h n 
(assumed pairwise different), subject to a primitive, reversible Markov process defined 
by a matrix M with left-P F- eigenvector p as described above. 

Then, the autocorrelation 7^ of oj exists with probabilistic certainty and has the 
form 7^ = Ylm&z u { m )^m with non-negative autocorrelation coefficients 

v[m) = (h|nM H |h) (21) 
form G Z. In particular, u(0) = (\h\ 2 ) = Ym=i Pi\^i\ 2 - 



PROOF: This is another application of Birkhoff 's pointwise ergodic theorem [f44| , Thm. 
2.3], this time applied to the case of an ergodic Markov system as described in pi , 
Sec. 1.2 D]. The setup is parallel to that of the Bernoulli system, only the invariant 
measure (defined via cylinder sets) is different, and this accounts for the different 
ensemble average. The latter is calculated as (for m > say) 

u(m) = V hi Pi Mi i Ma -...-Mi t \ = (h\UM m \h) . 

Due to (0), we also have 

(h\UM m \h) = (^|(M*) m n|/i) = (h\(UM m Y\h) = (h\UM m \h), 



which shows that v(m) = v(m), so u(m) is real. The case m < is analogous. But 
we also have 



v\m) 



(h | -(I1M H + (M*) H n) | h) > 0, (22) 
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because it represents a quadratic form with a non-negative real symmetric (hence 
Hermitian) matrix. Finally, m = gives the result for v(0). □ 

From the last Proposition, it is possible to derive the diffraction. Observe that, 
with c = Ph and r > 0, we have 

u(r) = (c\S r \c) = (c\S r \c) + (c\S[\c) , (23) 

and that So is a projector, i.e. Sq = S . With the explicit form of So given above, we 
thus obtain (r > 0) 



i>{r) 



2 

2 



+ (c\S[\c) = \(h)\' + (c\S[\c) . (24) 
With i/(0) = (\h\ 2 ) = | 2 + ((\h\ 2 ) - \ (h)\ 2 ), the autocorrelation gets the form 

oo 

7. = |W| 2 ^5 m + ((|h| 2 )-|(h)| 2 )5 + E( c l 5 il c )^ + M- (25) 

meZ r=l 

Before we continue, let us point out that the Bernoulli case is a special case of 
this, namely — 1/n, hence p 1 = . . . = p n = 1/n, S = M and Si = 0. In this 
limit, (|25| ) gives back the corresponding result of Proposition |I[ Also, one can treat 
more general cases of Markov chains, compare ]17[ vol. 1, Ch. XV] for details, and 
Markov chains of higher order (or depth). This becomes technically more involved, 
but we think that the essential flavour is obvious from the situation discussed here. It 
should also be clear how to make the other examples of the crystallographic literature, 
see pi] ], rigorous this way. Finally, let us remark that the similarity of the (binary) 
Markov chain to the ID Ising model is anything but accidental, and that the form of 
the autocorrelation is closely related to the solution of the Ising model by means of 
transfer matrices, see |58) and 0, Ch. 3.2]. 



The three terms on the right hand side of (|25|) can now easily be Fourier trans- 
formed. The first, by means of Poisson's summation formula (|i~3|), gives the Bragg 
part. The second results in a constant continuous background, (\h\ 2 ) — \ (h)\ 2 , as in 
our previous example. To calculate the Fourier transform of the third term, we can 
employ Neumann's series [46], p. 191] twice, because with Si also exp(±27r2/c)S , i is a 
contraction, for arbitrary k G KL This gives an absolutely continuous contribution 
which depends on the wave number k. Observing (\h\ 2 ) = (h\P 2 \h), one can combine 
the /c-dependent part with the first term of the constant part. This finally gives 
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Theorem 1 The diffraction spectrum of the Markov system of Proposition [I] exists 
with probabilistic certainty and is given by the formula 



% = \(h)\ 2 -^ + (%) ac , (26) 

where uj z = ^2 ke % o~k is the Dirac comb of the integer lattice. The absolutely continuous 
part (7^)00 is represented by the Z-periodic continuous (hence bounded) function 

^»-< ft ' S-2j(2/4s, + S ? P ' fc >-'W' 2 - < 27 > 
with obvious meaning of the quotient. □ 

Let us add that an explicit calculation of f(k) is easily done by means of the 
orthonormal eigenbasis of 5*. If Ph = ^™ =1 ^bj, one finds |/?i| 2 = \ (h)\ 2 and hence 



j=2 



which, for n = 2, coincides with the result of 0, [58fl . Note that the diffuse background 
corresponds to a positive entropy density which is given by s = — J2i Vi^ij l°g(^%)> 
53, Ch. 5.3, Ex. 3.5]. Let us also mention that (|27|) can also be rewritten as 



sec 



W = < h ' P l-2 1 c;K + ^ P ' h >- (29) 

from which f(k) > can easily be seen also in operator form. 

For a graphical illustration of the diffraction, we refer to f58fl . The important 
observation is that, in the presence of an interaction, the (now structured) diffuse 
background is attracted or repelled by the Bragg peaks according to the interaction 
being attractive or repulsive. This is a general qualitative feature of diffuse scattering, 
and will reappear in our later examples. 



ID random tilings 

Let us now change our point of view and consider a Bernoulli system not in the 
scattering strength h but rather in the distances. To keep things simple, consider the 
case of placing two intervals, of length u and v (both > 0), with probabilities p and 
q = 1 — p (hence, the entropy density per interval is s = — plog(p) — glog(g)), and 
assume that we have a unit point mass always at the left endpoint of each interval. 
Here, we have to distinguish the cases where a — u/v is rational or irrational. 
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Proposition 3 Consider the ensemble of binary random tilings of intervals of length 
u and v, with probabilities p and q = 1 — p, pq > 0. Let, for each such random tiling, 
A be the point set defined through the left endpoints of the intervals. Then, the natural 
density of A exists with probabilistic certainty and is given by d = (pu + qv)~ x . 

If io = oo A = XLeA $x denotes the corresponding stochastic Dirac comb, the auto- 
correlation 7^ of u also exists with probabilistic certainty. It is a pure point measure 
that is supported on the set 

A = {mu + nv \ m, n G Z and mn > 0} , (30) 

and, with z m ^ n := mu + nv, it is given by 

•& = <« E e'C.V^'-'w,^,^- pi) 

jv=_oo e=o v J 

For a = u/v irrational, it is of the form 7 W = YlzeA v(z)5 z , where z G A has 
a unique representation z = z m ^ n = mu + nv with m,n G Z and mn > 0. The 
corresponding autocorrelation coefficient is then given by 

v{Zm,n) = d( |n l + l m V m| S H - ( 3 2) 

V \ m \ J 

PROOF: Let A be a random point set according to the assumptions. If u = v, u A is 
the Dirac comb of a lattice and the statement is trivial. So, let us henceforth assume 
that u ^ v. 

If we view the random tiling ensemble as a Bernoulli system in the symbols u, v 
with attached probabilities p, q, each tiling is a sequence x = (x i ) ieZ with X{ G {u, v}. 
Since u, v also code the length of the intervals, the average distance between two 
consecutive points of the corresponding point set A is the limit of ^j-j- YH=- n x * as 
n — > oo which almost surely exists and, again by Birkhoff 's pointwise ergodic theorem, 
is given by (pu + qv). The density d is then clearly the inverse of this, as stated. 

The possible differences between two points of A are clearly given by A of Eq. (|30|), 
and this set is discrete and closed. To establish the existence of the autocorrelation, it 
is thus sufficient to show that its coefficients exist. Let z = mu + nv be in A. Although 
the representation of z need not be unique (e.g. if a is rational), there is no other 
representation with N = m + n intervals because we have excluded the case u = v. 
So, starting from an arbitrary point x G A, + 1 different points can be reached 
by adding A^ intervals (to the left or to the right, according to the sign of z), and 
the corresponding probabilities follow a binomial distribution, because the system is 
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Bernoulli. If we pick a single sequence, and determine the average frequency to reach 



the point x + z from x in N steps, the strong law of large numbers [17, vol. 1, Ch. X.l] 
tells us that this frequency, almost surely, converges to the corresponding probability, 
i.e. we find the limiting frequency 



n\ + \m\ 
\m\ 



p \m\ q \n\ 



Here, we have averaged over the number M of starting points x in a finite piece of 
the sequence A, and considered the limit M — > oo. The corresponding contribution to 
the autocorrelation coefficient, however, is defined as a volume- averaged limit. This 
gives a prefactor that is the average number of points of A per unit volume which is 
the density d. It exists with probability one as shown above. This establishes (|3~1~D . 
The full autocorrelation coefficient is now given by 

"M = *s.h £ 1 = d £ ( W M H ) pW9W ' (33) 

xeA, \x\<r m,neZ,mn>0 
x+zeA mu+nv=z 

which clearly also exists with probability one. This shows the existence of j w . 

Let us finally assume a ^ Q. Then, z m ^ n = z m i >n i implies m — ml and n = n', so 
that the only possibility to fill this distance is by m intervals of length u and n of 
length v, the remaining freedom just being the order in which this is done. So, the 
sum in the previous equation reduces to one term, the one given in (|3"2"|). □ 

Let us add two remarks. First, the autocorrelation could also be worked outQ by 
means of the renewal theorem |T7|, vol. 2, Ch. 11]. This would have the advantage 
of also being applicable to gases. However, for our case, one has to pay attention to 
convergence questions wherefore the derivation is not shorter. Second, the argument 
given for the existence of the average distance between two consecutive points of A can 
easily be modified to calculate that the average distance bridged by A" consecutive 
intervals is given by N(pu + qv), because the system is Bernoulli. This average, 
however, can now also be calculated as the weighted sum over the possibilities to fill 
A^ steps by m intervals of type u and N — m of type v , i.e. we obtain the identity 

( p ) P^'" ^ u + ( N ~ £ » = N(pu + qv) (34) 



which can also be checked explicitly by induction. It rests upon p + q = 1, and the 
binomial formula for (p + q) ■ It can also be understood from the first moment of the 
binomial distribution, calculated as derivative of its generating function. 
2 We thank A. Martin-Lof for pointing this out to us. 
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To understand the diffraction, we have to determine the Fourier transform of r y lJ . 
To do so, it is advantageous to write the tempered measure 7^ as a weak limit of 
tempered measures with compact support, i.e. to write 7^ = liniAr^oo n N where 

N \n\ ,v 

= d E E Y n m 1 >^ H " m w m , w _ m) (35) 

n=-Nm=0 V ' 

with z mn = mu + nv as before. It is evident that this sequence of measures con- 
verges weakly to 7^ of fl3~T|), and the support of fi N is certainly contained in the 
interval [— w,w] where w = N m&x(u,v). So, due to the Paley- Wiener theorem j|5|, 
Thm. IX. 12], the Fourier transform fi N is naturally represented by an entire analytic 
function, g^{k). Also, since the Fourier transform is continuous, the convergence of 
H N — > /i = 7^ implies that of fi N — > fi = 7 W . Note, however, that the n N are, in gen- 
eral, not measures of positive type, whence the jl N are not positive measures. They 
are signed measures though, as we shall see shortly, and one could decompose them 
as fi N = - fj,~ with fl% = l(\fi N \± £i N ). Then, jl% —> % and /% — >■ as N -> 00, 
but later calculations would be more complicated wherefore we prefer to work with 
the signed measures rather than with the positive measures fx^. 

Theorem 2 Under the assumptions of Proposition the diffraction spectrum con- 
sists, with probabilistic certainty, of a pure point (Bragg) part and an absolutely con- 
tinuous part, so 7^ = (tJpp + {l w )ac- If & — u/v, the pure point part is 

tfj = d 2 ■{ J * V ' (36) 

where, if a G Q, we se£ a = a/6 wift coprime a,b G Z and define £ = u/a = v/b. 
The absolutely continuous part (7 w )ac can 6e represented by the continuous function 

d ■ pqsm 2 (iik(u — v)) 
psm 2 (7iku) + g sin 2 (7rA;n) — pq sin 2 (irk(u — v)) 

which is well defined for k(u — v) G" Z. //; nas a smooth continuation to the excluded 
points. If a is irrational, this is g(k) = for k(u — v) G Z wift 7^ ana 1 

d-pq{u-v) 2 pq(u-v) 2 

0(0) = — 2-; — 2 ? « = rf 7 — , — \2- 38 

prr + qv — pq[u — u ) z (pu + gt>J z 

For a = a/b G Q as above, it is g(k) = for k(u — v) G Z ; £mi b^Z (or, equivalently, 
kv G" Zj, and 

*(*) = ^f^f (39) 
v y (pa + qb) 2 v y 



= — 2/ 7 x I ■ 27 , x • 2 /_,./_. ~ ' ( 37 ) 



16 



for the case that also ku G Z. 

PROOF: We employ the sequence of measures fi N introduced above which converges 
weakly to A direct calculation shows that the tempered measure jl N is represented 
by the analytic function 

N 

g N {k) = d (pe~ ssn{n)2mku + qe- sgn(n)2nikv ) H . (40) 

n=-N 

If we define r(k) = pe~ 2 " + qe~ 2 v , it is clear that this is a complex number inside 
the (closed) unit circle, i.e. r(k) = Re 1 ^ with < R < 1 and <fi G [0,2tt). This results 
in 

g N {k) = d (1 + 2 J2 Rmcos ( m( ^)] > ( 41 ) 

which shows that g N {k) represents a sequence of signed (or real), but not necessarily 
positive, measures. Their limit, however, is positive. 

Let us first check where the sequence of functions converges pointwise. We have 
r(—k) = r(k). Then, by the triangle inequality, \r(k)\ < p + q = 1. Also, we have 
|r(/c)| 2 = 1 — 4pqsm 2 (nk(u — v)), and thus |r(fc)| 2 = 1 if and only if p — 0, p — 1, or 
k(u — v) G Z. We have excluded the trivial cases p = and q = by our assumptions 
(they correspond to periodic chains, up to defects of density zero). So, if k(u — v) G" Z, 
the geometric series in fl40|) actually converges, with limit 

which immediately gives the expression in ( j37[) . In particular, the denominator is 
always different from for k(u — v) Z. 

It is not difficult to check that g(fc) has a continuation to points k with k(u — v) G Z. 
Consider first the case a irrational. Then, for k ^ 0, the denominator of (|37|) is never 
0, and = is the correct continuation. The case k = requires twice the 
application of de l'Hospital's rule, and gives the value of g(0) of (p8|). Next, let 
a = a/b with coprime a, b G Z. If k(u — v) G Z, then few G Z if and only if /ct> G Z. So, 
if few G" Z, we are back to the case where g(k) = is the correct continuation, while 
ku G Z, again with de l'Hospital's rule, gives the extension stated in (^). This in 
particular demonstrates that g(k), with the appropriate continuation, is a continuous 
function. It is also a positive function, as can easily be checked, and thus represents 
an absolutely continuous positive measure. 
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So, we have shown that g^{k) — g(k) tends pointwise to 0, as N — > oo, for all k 
with k(u — v) $ Z. The convergence is actually uniform on each compact interval that 
does not contain any of the exceptional points, as is clear from (f42|) . The latter form a 
ID lattice of spacing l/(u — v), and any singular part of ^ w must thus be concentrated 
to this set. Since the latter is uniformly discrete, the singular part cannot be singular 
continuous, but at most consist of point measures. 

We now have to check what happens with g^{k) for k{u — v) G Z, where the 
sequence of functions does not converge. First, let a be irrational, but k 7^ 0. Then, 
it is impossible to have r(k) = 1, because this would imply b£Z and hence teGZ 
- a contradiction to a G" Q. But then, with r(—k) = r(k) and \r(k)\ < 1, we get 



2\9 N (k)\ 



< 1 



< 1 



A? 



n=l 



1 + 2 



|r(A;)| |1 — r(k) 
|1 -r(k)\ 



N I 



N 



1 + \r(k)\ 
\l-r(k)\ 



< 1 



4 



1 + 



\l-r(k)\ 
2 



(43) 



a/p sin 2 (nku) + g sin 2 (7r/;;f ) — pq sin 2 (ixk(u — u)) 

which, for each fixed 7^ 0, has a denominator 7^ 0. So, the sequence g N (k), and 
hence fliOl) , stays bounded in this case, even though it does not converge. With 
the previous result on g, this means that the sequence g N is uniformly bounded on 
each closed interval that does not include 0. Hence, 7^ cannot be singular at the 
exceptional points k with k(u — v) G Z unless k = 0. The analogous argument applies 
if a = a/b G Q, as long as ku G" Z: the sequence g N is uniformly bounded on each 
closed interval that does not include any of the exceptional points with ku G Z, which 
we will call singular from now on. 

So, the remaining cases are k = for a irrational resp. ku G Z for a rational. For 
such singular k, ( ^D|) gives 



g N (k) = d(2N+l) 



(44) 



and this means that g N {k) diverges for these k, always with the same rate, and the 
divergence is proportional to the system size. To make this precise, consider first 
a = a/b with coprime a, b G Z. The singular points are then k = laju = Ibjv for 
£ G Z. Define L(N) = \_N(pa + qb)\, where [^J denotes the integer part of a positive 
x, and set 



h N (k) = d 2 i 



L{N) 
=-L(N) 



(45) 
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with £ = u/a = v/b. Clearly, h N (k) = d(2N + 1) + 0(1) for all singular k. The form 
of h N (k) is fixed by the requirement that both the height and the width of the finite 
approximations to the point measures at the singular points equals that of g N (k), up 
to lower order terms. One can now show (though we will skip the details here) that 
g N {k) — h N (k) is bounded on each compact interval even if it does contain singular 



values of k. On the other hand, we know from (|14D that 



hm h N (k) = d 2 £ e ' 2mkX = d " E S v ■ ( 46 ) 

This then proves that g N {k) converges to a point measure concentrated on |Z plus 
the ac measure derived above. Similarly, one deals with the case a irrational, but 
k = (e.g. by taking a suitable sequence of rational cases while letting £ — > 0). 

So, for all singular k, the structure factor converges to d 2 , i.e. (7 a; ) PP ({^}) = d 2 , 
and we obtain the result given in Eq. (|36|). 

Consequently, the positive measure 7^ has the decomposition claimed, and the 
absolutely continuous part can be represented by the continuous function g{k) (with 
the appropriate continuation to all k). This is the Radon- Nikodym density p6| , Thm. 
1.19] which is uniquely determined almost everywhere. □ 

Let us mention that the function g(k) appears simpler than it is — if one tries a 
selection of different parameters and produces some plots, one quickly realizes that 
it actually shows some "spiky" structure (though it is smooth), and the way it does 
depends rather critically on the nature of a = u/v. For example, if a = r (the golden 
ratio), a rather regular pattern emerges, and localized bell-shaped needles of increasing 
height appear at sequences of positions that scale with r. This is reminiscent of what 
happens in perfect Fibonacci model sets. If, however, a is transcendental (e.g. ir), 
much more pronounced needles appear, but at rather irregular positions. This is a 
clear consequence of how these numbers can be approximated by rationals, and an 
analogous phenomenon is well known in the approximation of irrational numbers by 



finite continued fractions [|TT|, Ch. VII]. 

We have discussed the case of a binary random tiling in detail, to make the structure 
as transparent as possible. It is clear that one can treat, with the same methods, also 
the case of a random tiling with n tiles of length (%,... ,u n ) = u, u { > 0, and 
attached frequencies (p 1( . . . ,p n ) = p, •p i > 0, Y^i=iPi = 1- us state the results in 
an informal way, as they are extremely parallel to what we discussed above. Viewing 
this system again as a Bernoulli system in the symbols u i reveals that the mean free 
path between two consecutive points of A is, almost surely, given by p ■ it, and the 
density is then d — l/(p ■ u). To simplify the following formulas, it is advantageous 
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to adopt standard multi-index notation. So, m = (m 1 , . . . ,m n ) is a vector of non- 
negative integers, \m\ 1 — m 1 + ■ • • + m n its 1-norm, and p m = p™ 1 ■ . . . ■ pn n . Also, 
we shall need the multinomial coefficient 

m J m 1 !m 2 ! • . . . • m n \ 
where N = \m\ 1 . 

Let A be the vertex set of a random tiling of this kind. The natural autocorrelation 
of cu = ou A exists with probabilistic certainty, and has the form 7^ = EzeA v i z )$z 
where A = {±z | z = m ■ u} and the autocorrelation coefficient is given by 



nrtrt.11 "y N / 



(48) 



This is the previous result with the binomial structure replaced by a multinomial one. 
In particular, we also get an analogue of Eq. (|34j), namely 



(jP m ('rn-u) = N{p-u). (49) 

iTni 1= jv y™ 1 / 

The autocorrelation measure 7^ can again be approximated by a weakly converging 
sequence of measures fi N , and their Fourier transform now reads 

N / n \ H 

Aiv = d- ^^.e-^ 2 ^ . (50) 

m=-N \j=l J 

The analysis of pointwise convergence then reveals once again that the diffraction 
spectrum consists of a pure point part and an absolutely continuous part. 

The absolutely continuous part of the diffraction, (7J00 is represented by the 
continuous Radon-Nikodym density 

g W = ^ ■ 2/ ^ V ^TTT ( 51 ) 

z2j Pj sm - z2j<e PjPe sm (nkiUj - u e )) 

which is well defined as long as not all • — u e ) are integer. If they are, the contin- 
uation is to g(k) = if ku U 1 and otherwise to 

, k , = d -J2 J< ePjPe(u j -u e ) 2 = d ■ J2j<i PjPt(Uj - u ef ^ 
9 Ej Pjtf ~ J2j<e PjPi( u j ~ u e) 2 {p-u) 2 
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The Bragg part of the diffraction is determined by the condition that, whenever 
ku G Z n , k is a pure point, and results in a contribution of d 2 5k to (xJpp- Thus, 
if u = £(a 1; . . . , a n ) with a { G Z and gcd(a 1 , . . . , a n ) = 1, the condition fctt G Z n is 
equivalent to A;£ G Z, and we get a ID lattice Dirac comb, 

(%)pp = d 2 ■ • (53) 

If, however, at least one quotient u i /uj is irrational, we only get the trivial point part, 
d 2 5o, while all other peaks are extinct. 

One reason for this rather detailed discussion will become apparent shortly when 
we use this to describe a class of very simple stochastic tilings in higher dimension. 



Intermezzo: Stochastic product tilings 

With the use of ID random tilings one can construct a particularly simple class of 
stochastic tilings in higher dimension, by simply taking one ID random tiling per 
Cartesian direction and considering the space filling by cuboids obtained that way. 
The prototiles are thus cuboids whose edges in direction j is any of the possible lengths 
of the jth ID random tiling used. Since the diffraction theory of these objects is 
essentially an exercise in direct products, but nevertheless quite useful and instructive, 
we describe the result in an informal way. Note, however, that the entropy density of 
these tilings is zero if D > 1, so that they are no ordinary random tilings in the sense 



of [g|] or g7|. 



Consider D different ID random tilings, and the corresponding point sets Aj, 1 < 
i < D, characterized by vectors of possible tile lengths ttW and frequency vectors pW. 
The total number of tiles in each case may be different, and is given by rij. Let us 
now consider the Cartesian product 

A = Aj x . . . x A D = . . . , x D ) | x t G Ai} , (54) 

which is the vertex set of a stochastic tiling in D dimensions whose prototiles are the 
n l - . . . -n D cuboids obtained as Cartesian products of the intervals vf? with 1 < i < D 
and 1 < j < Hi. The sets A are thus all of finite local complexity, and we have 

A = A- A = {z=(z 1 ,... ,z D ) \ Zi e (55) 

with probability one, where Aj = A,, — A,,. Also, the density of A exists again with 
probabilistic certainty, and is given by d = d 1 ■ . . . ■ d D where di = (pW ■ w^) -1 as 
derived in the previous Section. 
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The product structure of A also implies that u = J2 x ^a^x is a product measure 
(or, which is equivalent in this case, the tensor product of distributions, see |54|, Chap. 
IV]), i.e. we have 



D D 

-1 \xiSA 



ik' = n e<j? , (se) 



where 5^ is meant as a ID Dirac measure acting in the space along the ith coordinate. 
It also follows, with probabilistic certainty, that the autocorrelation 7^ exists and is a 
product measure, too. To prove the existence of the corresponding coefficients v{z), it 
is easiest to take averages over cubes rather then balls, i.e. to use the definition of [p6fl . 
Here, this gives the same limit as our definition of the natural autocorrelation. One 
obtains v{z) = Yl?=i ^H^) > where v®(zi) is the coefficient of the autocorrelation 
attached to A,-. So we have 



7« = TT Yl "®(*<)8$ ■ (57) 




Finally, let us consider the diffraction spectrum j u . Since it is the Fourier transform 
of a product measure, it is a product measure itself |54], Thm. XIV], and we thus obtain 

D D 

% = n%o) = n ((%w)pp+ {%(.*)) ac) > (58) 

8=1 1=1 

where the 7 (i ) are determined through Theorem ^. In particular, the absolutely 
continuous (pure point) part of 7^ is precisely the product of the ac (pp) parts of the 
7 w , while all other combinations result in singular continuous components — though 
the meaning of this will require some thought. The sc property can be seen from the 
fact that terms in the expanded product which contain at least one component of 
each kind are concentrated to a support of vanishing Lebesgue measure, but contain 
no pure points themselves — whence they must be singular continuous relative to 
Lebesgue measure. This also agrees with the common, intuitive scaling picture: a 
term with m ac-components and D — m pp-components would show intensities that 
stem from amplitudes (or Fourier-Bohr coefficients) which show a finite-size scaling 
with i/ n / 2 i J D - rn = [L D Y where L is the linear system extension and 

i </»=!-£<!. (59) 

Here, /3 — | and (3 = 1 correspond to the cases of ac andpp part, respectively, compare 
the discussion in [^6], Sec. 6] and f29fl . 
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However, this notion of singular continuity is to be taken with a grain of salt. 
All we have constructed here are product measures, and such objects would perhaps 
not qualify to be singular continuous in a "generic" sense. They do show up in liquid 
crystals though, compare the discussion on the nature of their long-range order in 



In particular, they appear in Danzer's aperiodic prototile in 3-space [13, Sec. 4], which 
can be seen as a toy model of a smectic C* liquid crystal. Nevertheless, these simple 
product tilings show that one has to expect a larger variety of spectral types in higher 
dimensions, and that already in planar cases the appearance of singular continuous 
components should be typical. Also, one can easily construct examples with all three 
spectral types present. More specific and genuine random tilings, however, might 
bypass this, as we shall see in the next Section, although that should not be considered 
generic. 

Two-dimensional random tilings 

Let us now move to planar systems, where we will mainly consider two illustrative 
examples, namely the classical random tilings consisting of dominoes and lozenges. 
Because of their symmetries, we call them crystallographic random tilings. Though 
still only two-dimensional, they are of practical relevance because of the existence of 
so-called T-phases (see || and references therein) which are irregular planar layers 
stacked periodically in the third direction due to a very anisotropic growth mechanism. 
It is thus appropriate to investigate the diffraction spectrum of a single layer obtaining 
then the complete spectrum once again as a product measure, compare the previous 
Section. 

Unfortunately, already the treatment of planar systems is a lot more involved than 
in the ID case. Although we will have to deal "only" with the action of Z 2 , we cannot 
directly apply standard results of ergodic theory as above, because we first have to 
establish the ergodicity of the measures involved. Even for the two simple systems 
we shall discuss below, the rigorous classification of invariant measures is only in its 
infancy, see |TI| and the discussion in [37]. Fortunately, the investigation of invariant 



equilibrium states, which form a subclass, is well developed [|TD|, |3D], [55] . If combined 
with certain results of statistical mechanics |5D[ , this allows for a determination of 
extremal states, which will be unique in our examples. They are ergodic and thus 
admit the application of Birkhoff's pointwise ergodic theorem in its version for Z 2 - 
action, see e.g. [^, Thm. 2.1.5]. Let us summarize the key features in a way adapted 
to our later examples. 
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Preliminaries 



A tiling u of a region^ A C M. d (with positive volume vol (A)) is a countable covering 
of A by tiles, i.e. by bounded closed sets homeomorphic to balls, having pairwise 
disjoint interiors and non-vanishing overlap with the region A. In our case, the tiles 
are translates of finitely many prototiles. We deal with free boundary conditions in 
the sense that the tiles may protrude beyond the boundary. Thus, the boundary of A 
does not impose any restrictions of the kind known from fixed boundary conditions or 



exact fillings of given patches, compare [3q. Two coverings of the same region A are 



called equivalent if they are translates of one another. For further conceptual details, 
we refer to ||7| . 



The examples discussed below belong to the class of polyomino tilings, where the 
prototiles are combinations of several elementary cells of a given periodic graph G. 
These tilings can be described as polymer models |[47|| , and as dimer models in our 
case. 



Figure 1: Representation of a random tiling on the dual cell complex. 



A dimer is a diatomic molecule occupying two connected sites of a graph. A graph 
is close-packed if all sites are occupied precisely once. In Figure [l], the one-to-one 
correspondence between the tiling on a (periodic) graph and its close-packed dimer 
configuration on the dual cell complex (the so-called Delone complex) is illustrated 
for the domino tiling. The scatterers (Dirac unit measures) are placed in the centre 
of the tiles, resp. dimers. 

In the space Q of all tilings, two elements are close if they agree on a large neigh- 
bourhood of the origin. Q is compact in this topology. The group Z 2 of translations 
acts continuously on Q (in an appropriate parametrization) because of the periodicity 
of the underlying graph. 

We use a grand-canonical setup where we assign equal (zero) interaction energy and 
a finite chemical potential fa or activity Zi = e Mi to each of the M different prototiles 

3 We tacitly assume that any such region is sufficiently nice, i.e. it should be compact, measurable 
and simply connected. 
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(setting the inverse temperature (3 = 1). For fixed prototile numbers n%, . . . ,Um, 
let us denote the number of nonequivalent A-patches that use prototiles of type i 
by g/i(ni, . . . ,Um)- The grand-canonical partition function is given by the following 
configuration generating function (/i = (/i l5 . . . , [X M )) 

Z M = Yl 9 ^ nu ■■■,n M )z^ 1 ■ z^ M . (60) 

n v ...,n M 

To adapt the usual language of statistical mechanics, let to G Q\ be a tiling of a finite 
region A which is now positioned relative to a fixed lattice, Z 2 say. An interaction is 
the translation invariant assignment of a continuous function $(o>) G to every uj, 
so that a shift of A by a G Z 2 results in $(u; + a) = (T a $)(co>), where r a denotes the 
corresponding shift (cf. |50], |55|). To avoid repeated counting of contributions to the 
interaction energy, $ represents only the basic interaction of tiles in uj that are not 
already included in the interaction of subsystems. For uj C Q A , we now define the 
Hamiltonian Hf(u>) by 

H*( U ) = $ ^') ' ( 61 ) 

U) 1 < UJ 

where the sum runs over all sub-patches of to. This is well defined due to the re- 
strictions on $ mentioned before. Let B be the Banach space (with norm ||.||oo) of 
interactions $ subject to the restriction 

|||*||| := ^)l!*ffi < co, (62) 



where the sum is over all tilings covering the origin and \u\ denotes the number of 
tiles (cf. also 0, App. B]). Let us now, for simplicity, assume that each tiling to of A 
has the same total number of tiles, \uo\ = N A . The pressure (negative grand-canonical 
potential) per tile is then given by 

Pa(*) = aT E e ^* H = ^log^ A ($). (63) 
A ujen A A 

If the extra assumption is not fulfilled, A^ A is the average number of tiles. This is a 
reasonable definition as long as we take the limit A — ► oo in the sense of van Hove, 
see |)D] for details on this concept. But then, for $ G £>, the thermodynamic limit 
p($) = liniA^ooPA^) exists and is a convex function of $, compare PD|, The 
interaction in our models is simply the self-energy (— //j for prototiles of type i) and 
thus included in B. 
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A state v of the infinite system is a Borel probability measure on Vt. A state v is 
called invariant or translation invariant if it is invariant under the action of J?. The 
set of all invariant measures, A4 1 , is compact and forms a simplex. 

Following Ruelle's presentation 0, Ch. 7.3], let D be the set of all $ G B such 
that the graph of p has a unique tangent plane at the point ($,p($)). If $ G £>, there 
exists a unique linear functional a* in the dual B* of i3 such that 

p($ + ^) > p($)-a*(tf). (64) 

The mean densities pf(p) = ^2i(A) > * = lj--->-^> °f the different prototiles in the 
ensemble can be computed as functional derivative of p\(p) with respect to p (( • )a 
denotes the finite-size average for given chemical potentials pi, . . . ,Pm)- For p G B 
and A — > oo in the sense of van Hove, we haveQ 

M 

lim ^Tpf^ = a"(/2) (65) 

i=l 

for all p G B, and 

A dp(p!, . . . , p M ) (aa , 

hm p. = ft = tt , (66) 

a^oo a pi 

with Yji=iPi = 1- Since the chemical potentials and also the conjugate (mean) densi- 
ties do not form an independent set of macroscopical parameters, we may choose an 
independent subset p 1 , . . . , p k by setting p k+1 = ■ ■ ■ = p M = 0. This normalization of 
p leaves the densities invariant. 

The entropy per tile of a finite region tiling is defined as 



T7- u Mlogv A (u) (67) 



= N A 

where u A is the restriction of v to Sl A . For translation invariant measures and p G B, 
the infinite volume limit exists, giving u A — > v for A — > oo taken in an appropriate 



way, and the functional s is afline upper semicontinuous |j0[ Ch. 7.2]. According to 



Gibbs' variational principle |50|, Ch. 7.4], the pressure can then be calculated as 

p(p) = sup [s(u) - v{p)\ . (68) 
veM 1 

The measure for which this supremum is attained is called equilibrium measure. We 
may formulate the weak Gibbs phase rule [SO, Ch. 7.5]. 



4 Usually one defines A$>{uj) = XL'<w ^\^'\ e ®* ■ Restricting ourselves to the tiling models where 
the chemical potentials are the only interactions, we may identify and /i. 
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Theorem 3 Let D C B defined as above. 



1. If fi G D, the function v t— > s(v) — z/(/i) reaches its maximum p{ii) at exactly 
one point G M. 1 . 



2. If ll G D and a M G B* is defined by (p^l), then, for all fx G B, 

v*(fi) = of(fi (69) 

so that is the infinite volume equilibrium state corresponding to the chemical 
potential \i. 

3. J//iG D, is a l?-ergodic state and may thus be interpreted as pure thermo- 
dynamic phase. □ 

In what follows, we calculate the ensemble average of the correlations. Since the 
diffraction image is taken from a single member of the ensemble, the above theorem 
ensures that the typical member is self-averaging as long as the pressure is differen- 
tiable (no first order phase transition). 

Calculating the diffraction of our models consists essentially in calculating the 
corresponding dimers autocorrelation which we will base upon previous work of Fisher 
and Stephenson |]TB[ and of Kenyon p7[ . 



Kasteleyn has shown that for any finite planar graph with even number of sites, 
and also for any periodic graph with a fundamental cell of an even number of sites, a 
Pfafnanf] can be constructed which is equal to the dimer generating function. For this, 
one has to orientate the graph in such a way that every configuration is counted with 
the correct sign. In addition, every bond is weighted with the corresponding dimer 
activity Zj. The configuration function is then given by the Pfaffian of the activity- 
weighted adjacency matrix A. Although Kasteleyn's proof applies to arbitrary graphs, 
the calculations simplify considerably when restricted to periodic simply connected 
graphs as in our case. If we define occupation variables for a bond between sites k 
and k' 

(1, bond (k,k') occupied, 
(70) 
0, otherwise, 



we can state ITS 



3 A Pfaffian is basically the square root of the determinant of an even antisymmetric matrix, see 



App. E] or 1 42, Ch. IV.2] for an introduction. 
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Proposition 4 Let G be an infinite simply connected periodic graph with close-packed 
dimer configuration where each dimer orientation has density p i > 0. Let A be the 
invertible weighted adjacency matrix. If the dimer autocorrelation (joint occupation 
probability) exists, it is given by 

PaP = ( 71 ) 
= (Vh y )(Vu y ) - A k y A k k , {AT 1 . At},, - At},, A,}, ) . 

Proof: Let fj kh , = 1 — r] kk ,. Obviously, 

(Vk.KjlkpUf) = (V kak ' a )(v kpk ' l3 ) + (Vk a k' a Vk k' g ) - (Vk a k' a )(Vk f3 k' p ) ■ (72) 

The first term on the RHS of ((72|) depends only on the densities p a resp. p^ of the 
dimers that can occupy the bond (k a k a ) resp. {k^k'^). In the case of k a and k' a being 
connected by a bond of type i and only one bond of this type leading to each site, this 
would result in (r] k k /) = Pi, because we normalize with respect to the total number 
of dimers (and not to the number of sites or bonds as in [JTSj] ) . It remains to prove 



so we 



the equivalence of the second terms of (|7ll) and (72). This was shown in |18 
will just give an outline here. 

Here, {fj kak ' a ) is the (weighted) sum of all dimer configurations where the bond 
(k a k' a ) is not occupied, divided by the total (weighted) sum of configurations Z. 
Since Z = Pf(A), we define Pf(A) = Pi(A + E) as the "perturbed" Pfaffian counting 
precisely all configurations where (k a k' a ) is not occupied, i.e. all elements of E are 
zero except E k k > = —E k > k = —A k k > . Note that A, A and E are skew-symmetric 
matrices. So (fjk a k' a ) = Pi(A)/Z = Pf(I + A~ 1 E) = Pf(E) P^-ET 1 + A" 1 ), where 
the last equality holds only if E is invertible, I denotes the unit matrix of appropriate 
dimension. The same applies to (f]k a k' a Vk l3 k' ) but with four nonvanishing elements of 
E. The evaluation of the Pfaffian reduces to the calculation of a small determinant 
and yields the desired result. □ 

The calculation of Z and A -1 can be simplified considerably by imposing periodic 
boundary conditions. Usually, for the partition function for a toroidal graph, one 
needs four determinants differing from that with free boundary only in exactly these 



boundary elements ]53]. But in the infinite volume limit, these modifications do not 
change the value of a determinant (given by the product of the eigenvalues) as can be 
derived from the following result of Ledermann 



Lemma 1 If in a Hermitian matrix the elements of r rows and their corresponding 
columns are modified in any way whatever, provided the matrix remains Hermitian, 
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then the number of eigenvalues that lie in any given interval cannot increase or de- 
crease by more than 2r. □ 



If N = run is the number of sites or elementary cells, the number of eigenvalues per 
unit length, which is O(N), changes only by 0(\/ r N), which is negligible in the limit 
as iV — > oo (for details see also JJ3J). The same argument holds for A -1 . 

How can we calculate the elements of A -1 ? Since A is the adjacency matrix 
of a graph made up of a periodic array of elementary cells with toroidal bound- 
ary conditions and is therefore cyclic, it can be reduced to the diagonal form A = 
diag{Aj} by a Fourier-type similarity transformation with matrix elements S kk > = 
{mn)~ 1 / 2 exp(27ii(k 1 k[/m + k 2 k' 2 /n)). A -1 is now determined by 

(m,n) 

A k l = {SAT*S-i) M = ]T S kj X/Si,.. (73) 

j=(i.i) 

In the infinite volume limit, the sums approach integrals (Weyl's Lemma), and by 
introducing (p 1 = 2iTij l /m etc. we obtain 

2tt p2tt 



All, = i^/ j ^" 1 (ft.Wi) e " <( " l( * 1 "' lH *' W " y) ^ ( 74 ) 
Let us illustrate this by two examples, see p5) for another case. 



Domino tiling 

A domino is a 2 by 1 or 1 by 2 rectangle, whose vertices have integer coordinates in 
the plane. A tiling of the plane with dominoes is equivalent to a close-packed dimer 
configuration on the square lattice Z 2 . This model exhibits no phase transition. We 
assume finite, positive activities in order to have non-vanishing tile densities. The 
degenerate case will be treated separately. Labelling the sites by k = (A; 1 ,/c 2 ), and 
adopting from the various equivalent possibilities (see |TB[ for details) the technically 
most convenient choice of complex weights, one obtains the weighted adjacency matrix 

A(ki, k 2 ; k\ + 1, k 2 ) = —A{h\ + 1, k 2 ; fa, k 2 ) = zi, 
A(fa, k 2 ] fa, k 2 + 1) = —A(fa, k 2 + 1; fa, fa) = iz 2 , (75) 
A(fa, fa; k[, k' 2 ) = A kk i = 0, otherwise, 

with eigenvalues 

X((p 1 ,ip 2 ) = 2i (z x sin <p x + iz 2 sin cp 2 ) . (76) 
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After inserting in ( [7ip and separating real and imaginary parts, one gets (k' — k = 
r = ( Xl y)eZ 2 ) 

= ^J J ^sinV 1+ 4smV 2 ^ 2 (7?) 



with 



'o ^sm ^l + ^2 Sm ^2 

0, (x, ?/ same parity) 

M(x,y\tp 1 ,tp 2 ) = ^ — z x siny9 1 sinx^ cosy(p 2 , (x odd, ?/ even) (78) 

-iz 2 sin <^ 2 sin ?/<y? 2 cos xcp x , (y odd, x even) 

(compare with fl8|]). We introduce the abbreviation for the coupling function [[TjJ 

[x,y] = A- x {k x MXX), (79) 

obeying [x,y] = —[—x,—y]. We place the scatterers in the centres of the tiles or 
equivalently of the dimers. With (|71|) and (|72|), the joint occupation probability of 
two horizontal dimers with scatterers at distance r in the centre of an infinite lattice 
is given by 

P n (r) = | + c n (r) = & - z\ ([x, y} 2 - [x - 1, y}[x + 1, y}) , (80) 

(similarly for P 22 ) where only c 11 (r) (resp. c 22 (r)) depends on the distance r. For a 
pair of mutually perpendicular dimers, the possible distance vectors of the scatterers 
r + a, a = (—1/2, 1/2)*, are odd half-integer. The joint occupation probability is 

P 12 (r + a) = P 2 i(r-a) 

= El?*-i ZlZ2 ([ X ,y][ X -l,y+l]-[ X ,y + i\[ X -l,y]) . (81) 

The non-constant part is c 12 (r) = c 21 (r). Note that either the first or second part of 
the term in brackets vanishes because of the parity of x and y. The full autocorrelation 
for the positions of the scatterers is thus given by 

7. = E ( Pl1 ( r ) + P22 &r + * E ( Pl2 (»" + «) + ^21 (r + a)) 5 r . (82) 

rel? rez 2 
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Theorem 4 Under the above assumptions, with pip 2 > 0, the diffraction spectrum of 
the domino tiling exists with probabilistic certainty and consists of a pure point and 
an absolutely continuous part, i.e. 7^ = (tJ^ + ("y u ) ac , with 

(%)„ = \ E (px + (-l) h+k P2) 2 S (h>k) . (83) 

(h,k)& 2 

In particular, there is no singular continuous part. Furthermore, j u is periodic with 
lattice of periods {(h, k) G Z 2 | h + k even}. 

PROOF: The point spectrum can be calculated directly by taking the Fourier trans- 
form of the constant part of fl82"l) . This requires Poisson's summation formula and the 
convolution theorem, which leads to the phase factor (— V) h+k and to the periodicity 
claimed. The non-constant part of 7^ is determined by c u , c 12 and c 22 . To establish 
our claim, we will show that the formal Fourier series ^ rgZ 2 c i j(r)e~ 2mk r actually 
converge to /^-functions and thus represent absolutely continuous measures. 

The real coefficients c i; j(r) are essentially products of the form [x, y] 2 . One in- 
tegration in ( |77| ) may be performed explicitly. With standard asymptotic methods 
involving the Laplace transform, see the Appendix of [ BO | , one obtains the asymptotic 
behaviour in the limit of large x and y as 



1 



[x odd, y even), 



tt (z, 2 x)' 2 + (z 1 y) 2 ' 



— (y odd, x even) 



Thus, we obtain < Cy(r) = O ^ x2+ y 2+1 ^ < 00, since x and y have different parity. 
Note that the implied constant still depends on z 1 and z 2 - Now, e.g. by referring to 
Eqs. (6.1.126) and (6.1.32) of 0, one can see that J2^=o ^27=o r^+^+Tp converges 
(it actually is even less than 2 in value). Consequently, by Cauchy's double series 
theorem, J^rez 2 ( c *i( r )) 2 converges absolutely. So, the Cij(r) can be seen as functions 
in £ 2 (Z 2 ) and, by the Riesz-Fischer Theorem [^4[ Thm. 23.3], each of the Fourier series 
^2 re z 2 c i j(r)e~ 2mkr converges to a function in L 2 (R 2 /Z 2 ) in the L 2 -norm. The limit 
is independent of the order of summation, and convergence is actually also pointwise, 
almost everywhere. Holder's inequality gives L 2 (M 2 /Z 2 ) C L 1 (IR 2 /Z 2 ), and combining 
these periodic functions with the appropriate phase shifts as implied by ( |8"2"D results 
in a function with lattice of periods r = {(h, k) G Z 2 | h + k even} which is certainly 
in L 1 (IR 2 /r), so the Radon-Nikodym theorem ^| leads to the result stated. □ 
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REMARK: One may assign complex weights h\ y hi to the scatterers on the two 
dominoes without changing the spectral type. This results in the pure point part 



(%)„ = \ E \pih + (-i) h+k P2h 2 \ 2 8 ih , k) 

(h,k)ez 2 

This applies analogously to the next example. 
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Figure 2: Typical tiling for p 2 = 0.3 (left) and its diffraction image (right). The 
scatterers are located at the centre of the tiles. 

One has to be aware that the pure point part ( jB2|) does not display the correct 
(statistical) symmetry of the system. Away from the point of maximum entropy 
(which is p 1 = p 2 = 1/2), the tiling is no longer fourfold symmetric, as still indicated by 
the point part, but the twofold symmetry is only displayed in the diffuse background. 
This can be seen in Figure ^|. The pp part is calculated from the exact expression 
and the Bragg peaks are represented by white circles with area proportional to the 
intensity. The ac part was calculated numerically by means of standard FFT, because 
this is simpler than using the exact expression for the correlation functions. 

In the special case of only one domino orientation remaining, the scatterers distri- 
bution is a Dirac comb on a rectangular lattice. Using Poisson's summation formula, 
the diffraction spectrum is a Dirac comb on the reciprocal rectangular lattice. In 
such a limit, the diffuse background accumulates at the extra positions and converges 
vaguely to a point measure that completes the square lattice arrangement to the 
proper rectangular one. 
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Lozenge tiling 

A lozenge is a rhombus with side 1, smaller angle tt/3, and vertices in the triangular 
lattice T = A 2 /V2 with minimal distance 1. This tiling can be mapped on a dimer 
configuration on the honeycomb packing. The different tile densities are nonvanishing 



(a) (b) 



Figure 3: (a) Lozenge tiling and the dimer configuration of the honeycomb packing, 
(b) Elementary cell for the weighted adjacency matrix. 

if z i > \Zj — z k \, i,j, k pairwise different P7[| ; at equality, the system undergoes a phase 



transition of Kasteleyn type |35[ with only one lozenge orientation remaining. This 
trivial case shall be excluded in the sequel. With Dirac unit measures on the tile 
centres, the support of the scatterers is given by a Kagome grid of minimal vertex 
distance 1/2. The adjacency matrix then has entries that are 2 by 2 matrices them- 
selves, describing the elementary cells of the packing (see Figure |3|). By wrapping the 
graph on a torus, we may transform it to block diagonal form with elements 

= ( ■ ■ -(v^"+v-*+*)y (86) 

W 1 '™ \ Zl e 1 ^ + z 2 e 1 ^ + z 3 J 

If we introduce a coordinate system with x = (1, 0)* and y = 1/2(1, \/3)*, we can use 
the above notation for the difference vectors of the elementary cells. Denoting the left 
and right site of the elementary cell by L and R we see that [x, y) LL = [x, y] RR = 0. 
In the infinite size limit, the remaining matrix elements are given by 

[x,y\^^z 3 ] LR = ^ Zie -^ + Z2e -i, 2 + ^ 



(compare with |37||). Let v = e liPl and w = e ltp2 . Then 

r i 1 f v~ x w~ y dv dw , 00 . 

[x,y z v z 2 ,z 3 LR = — / ■ ■ — — . 88 

47T^ Js^-xS 1 z l v + Z 2 W + Z 3 W %W 
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As was already observed by Kenyon |37J for the isotropic case, the coupling func- 
tion has all the symmetries of the graph: interchanging v and w or the substitution 
(v,w) — > (w -1 ,™ -1 ) and combinations of these let the integral invariant, i.e. 

[X,V | Z 1 ,Z 2 ,Zg\ L r = [x, -x - y - 1 | z x , z 3 , z 2 ] LR 
= [y,x | z 2 ,z x ,z 3 \ LR = [y,-x-y -1| z 2 ,z 3 ,z x \ LR (89) 
= [— — x/ — 1, a; | z^z^z^lr = [-x - y - l,y \ z 3 , z 2 , z^lr . 



We evaluate one integration in ( |88"D explicitly for x < —1. The other values can 
be obtained by (|59D. This is a direct generalization of |S7J to the case of arbitrary 
activities. One gets 

[x.yW.z^z^LR = —{-z x ) x / uT* -1 (z 2 + z^w)^" 1 dw, (90) 

Ztt J e i Vo 

2 2 2 

with ip Q = arccos Zl 2 * 2 Zs . One easily finds the possible distance vectors of the scat- 
terers. Away from the phase transition points, we have for the constant part of the 
autocorrelation measure for the scatterers (one per lozenge) 

(iJconst = ((pI+ Pl+ Pl) 6 (x,y) 

V (*,»)er (91) 



x+l 2y + l\ 2y+l ' 

The reciprocal lattice T* is spanned by the vectors and ^0, -j=^ . Using 

(ffH) we can state 

Theorem 5 Under the above assumptions, the diffraction spectrum of the lozenge 
tiling exists with probabilistic certainty and consists of a pure point and an absolutely 
continuous part, i.e. % = {%) pp + {%) ac , with 

{l u )„=\ E ((-^V + t-l^ + PsfW)" ( 92 ) 

(h,k)eT* 

There is no singular continuous part, and 7 W is periodic with lattice 2T* . 



Proof: The pure point part is simply the Fourier transform of (|9*ID, again calculated 
by means of Poisson's summation formula. 

As before, we will now show that the remaining part of 7 W converges to a peri- 
odic L 1 -function and thus represents an absolutely continuous measure. Let us start 
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Figure 4: Typical tiling for p 2 = 0.24 and p l = p 3 = 0.38 (left) and its diffraction 
image (right). The scatterers are located at the centres of the tiles. 

with the case x < — 1 (and y arbitrary, but fixed) and show that [x,y\z 1 , z 2 , z 3 ]lr = 
O (|x| _1 ) as x — > — oo. The necessary extension to the complete asymptotic behaviour 
will later follow from Eq. flS9|). For now, and for fixed values of the activities in the 
admitted range, we have 

1 = I [^vIz^z^ZsIlrI < tt \ z 2 + z z w\~ x ~ 1 \dw\ 

[z\ + z\ + 2z 2 z 3 cos ■&) d-d. (93) 

We first bound /($) = (z 2 + z\ + 2z 2 z 3 cos "&)/zf by a straight line /($), i.e. we look 
for /($) < /($) on the interval [y? ,7r]. Note that /($) can only have one flex point 
at d = tt/2 in [0, it}. One has to distinguish two cases: 

1. cosv9 > (tt - <p„)siny? _- 1: Choose = (# - ip Q )f'(<p s ) + f(jp Q ). The 
condition on <p implies /i(tt) > /(7r). Differentiating /(#) — with respect 
to $ yields /'(i?) — = /'($) — f'ifo) an d this vanishes if $ = <^ (maximum) 
or if i? = tt — <£> (minimum). Thus /(i?) < 

2. cos<^ < (7r — y? ) siny? — 1: Choose ^(i?) = ^j > ~ 1 '$ + n ~^°^ n ^ (the line con- 
necting /(<p ) = 1 and /(7r)). Obviously J^l 71 ") — f'i 71 )- Because of the angle 
condition, we further have f((po) < /^(^o)- As there is only one flex point, we 
conclude that /(#) < / 2 (tf). 
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Choose /($) as or ^(i?) according to the previous distinction. If c = we 
get < c < (l - (22 ~2 3)2 ) < ^z^- Consequently, we can estimate 



-se-1 



TT^I < / (l-c(0-p o ))^~tft? 
Jip 

< - I (1 -t)^dt (94) 



c 







cfl — x) \\x 



Now, we can use the symmetry relations of (p9|) and obtain the asymptotic be- 
haviour [x, y\z 1 , z 2 , z 3 ]lr = O ((\x\ + whenever |r| — > oo, compare |37[| . The 



rest of the argument is very similar to the domino case and need not be repeated here, 
the periodicity statement follows again constructively. □ 

Let us remark that an analogous scenario, with the same type of result, occurs for 
the more complicated dart-rhombus random tiling, see [EH] for details. 



Addendum: The two-dimensional Ising model 

For the sake of completeness, we add an application of the probably best analyzed 
model in statistical physics, the 2D Ising model without external field. It may be 
regarded as a lattice gas on Z 2 , compare f55f , with scatterers of strength s^* G {1, 0}. 



The partition function in the spin-formulation (07^) € — 1}) reads as follows 

Z = J^exp J2%d) { K i a (i+i,j) + K 2*(ij+i)) > ( 95 ) 

M \(i,3) J 

where we sum over all configurations {er}. We consider the ferromagnetic case with 
coupling constants Ki = Jj/(fcsT) > 0, temperature T and Boltzmann's constant 
kg. The model undergoes a phase transition at k := (sinh(2i^ 1 ) sinh(2fT 2 )) _1 = 1. 
It is common knowledge that in the regime with coupling constants smaller than the 
critical ones (corresponding to T > T c ) the ergodic equilibrium state with vanishing 
magnetization m is unique, whereas above (T < T c ) there exist two extremal equilib- 
rium states, which are thus ergodic [53, Ch. III. 5]. In this case, we assume to be in 



the extremal state with positive magnetization m = (1 — k 2 ) 1 ' 8 . 

The diffraction properties of the Ising model can be extracted from the known 
asymptotic behaviour |42|, |59| of the autocorrelation coefficients. We first state the 
result for the isotropic case (K\ = K 2 = K) and comment on the general case later. 
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Proposition 5 Away from the critical point, the diffraction spectrum of the Ising 
lattice gas almost surely exists, is Z 2 -periodic and consists of a pure point and an 
absolutely continuous part with continuous density. The pure point part reads 

1. T >T C : (%)„ = \ J2kez? S k 

2. T < T c : (%) pp = p 2 £ fc6Z2 S k , 

where the density p is the ensemble average of the number of scatterers per unit volume. 

PROOF: First, note that s~ = (ov^-a + l)/2 and thus (^uj\) = m = 2p — 1, so p 
varies between 1 and 1/2. The asymptotic correlation function of two spins at distance 
R = ^/x 2 + y 2 (as R -> oo) is g| 



ci^, T>T C 



(^(o,o)^)) - {"\^ \-2 R/c2 ' ' ^ (96) 
vnf 1 + c 3 R2 , T < T c , 



with constants ci,C2 and C3 depending only on K and T, see also |39|, p. 51] and 
references given there for a summary. The pure point part {^J) pp results directly from 
the Fourier transform of the constant part of 7^ as derived from the asymptotics of 
( s (o,o) s (x,y)) = ((^o.o)^)) + 2m + l)/4. 

Here, already Y^( x y)& 2 e~ R ^ C2 / \^R and Yl( x y )^i? e~ 2R ^ C2 / R 2 converge absolutely, so 
we can view the corresponding correlation coefficients as functions in L 1 (Z 2 ). Their 
Fourier transforms (which are uniformly converging Fourier series) are continuous 
functions on M 2 /Z 2 , see f48| , §1.2.3], which are then also in L 1 (M 2 /Z 2 ). Applying the 
Radon-Nikodym theorem finishes the proof. □ 

Remark: At the critical point, the correlation function ( cr (o o)^^ y)) ^ s asymptot- 
ically proportional to i?~ 1//4 as R — > 00 [|59| , |39fl . Again, taking out first the constant 



part of 7 W , we get the same pure point part as in Prop. |5| for T > T c . However, for the 
remaining part of j u , both our previous arguments fail. Nevertheless, using a theorem 
of Hardy |], p. 97], we can show that the corresponding Fourier series still converges 
for k (jL Z 2 (a natural order of summation is given by shells of increasing radius). 

For k G Z 2 , where the Bragg peaks reside, the series diverges. But this can neither 
result in further contributions to the Bragg peaks (the constant part of 7^ had already 
been taken care of) nor in singular continuous contributions (because the points of 
divergence form a uniformly discrete set). So, even though the series diverges for 
k G Z 2 , it still represents (we know that 7^ exists) a function in L 1 (M 2 /Z 2 ) and 
hence the Radon-Nikodym density of an absolutely continuous background. On the 
diffraction image, we thus can see, for any temperature, Bragg peaks on the square 
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lattice and a Z 2 -periodic, absolutely continuous background concentrated around the 
peaks (the interaction is attractive). At the critical point, the intensity of the diffuse 
scattering diverges when approaching the lattice positions of the Bragg peaks. 

The same arguments hold in the anisotropic case, where the asymptotics still con- 



forms to Eq. (96) and the above, if R = R(x, y) is replaced by the formula given in [59 
Eq. 2.6]. The pure point part is again that of Prop. [5] with fourfold symmetry, while 
(as in the case of the domino tiling) the continuous background breaks this symmetry 
if K\ ^ K2. Let us finally remark that a different choice of the scattering strengths 
(i.e. ±1 rather than 1 and 0) would result in the extinction of the Bragg peaks in the 
disordered phase (T > T c ), but no choice does so in the ordered phase (T < T c ). 

Outlook 

The diffraction of crystallographic and of perfect quasi-crystallographic structures is 
well understood. The main aim of this article was to begin to counterbalance this into 
the direction of certain stochastic arrangements of scatterers, and to random tiling 
arrangements in particular. This requires a careful investigation of the diffuse back- 
ground and clear concepts about absolutely versus singular continuous contributions. 
Such problems are naturally studied via spectral properties of unbounded complex 
measures which we did for a number of simple, but relevant examples. 

While stochastic cuboid tilings would generically display a singular continuous 
contribution in the diffraction image, our results on the domino and the lozenge 
tilings show that no such contributions exist there. However, we do not think that 
this is a robust result. In fact, the natural next step would be an extension to planar 
random tilings with quasi-crystallographic symmetries such as 8-, 10- or 12-fold. Then, 
according to the folklore results, one should expect a purely continuous diffraction 
spectrum (except for the trivial Bragg peak at k = which merely reflects the existing 
natural density of the scatterers). This spectrum should then split into an ac and an 
sc part. 

But is the replacement of Bragg peaks by sc peaks significant? Scaling arguments 
and numerical calculations |3~3J indicate that the exponents of important sc peaks 



can be extremely close to 1 which means that their distinction from Bragg peaks is 
almost impossible in practice. Is it possible that this is one reason why structure 
refinement is so difficult for real decagonal quasicrystals? This certainly demands 
further thought, but it is not clear at the moment to what extent a rigorous treatment 
is possible. 

Finally, except for Bernoulli type systems 0] and for extensions by means of prod- 
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ucts of measures, we have not touched the "real" diffraction issues in 3-space. One 
reason is that we presently do not know of any interesting model that can be solved 
exactly (let alone rigorously), another is that, for random tiling models relevant to real 
quasicrystals, one expects bounded fluctuations (again, due to heuristic scaling argu- 
ments, see [^] and references therein). Consequently, the diffraction spectra should 
typically show Bragg peaks plus an absolutely continuous background. We hope to 
report on some progress soon. 
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